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Abstract 

We study Schwinger mechanism for gluon pair production in the presence of arbitrary time- 
dependent chromo-electric background field E a {t) with arbitrary color index a=l,2,...8 in SU(3) 
by directly evaluating the path integral. We obtain an exact expression for the probability of 
non-perturbative gluon pair production per unit time per unit volume and per unit transverse 
momentum d i^fi pT from arbitrary E a (t). We show that the tadpole (or single gluon) effective 
action does not contribute to the non-perturbative gluon pair production rate d i^p ^ . We find 
that the exact result for non-perturbative gluon pair production is independent of all the time 
derivatives d dt J^ where n = 1,2, ....oo and has the same functional dependence on two casimir 
invariants [E a (t)E a (t)] and [d a b c E a (t)E b (t)E c (t)] 2 as the constant chromo-electric field E a result 
with the replacement: E a — > E a (t). This result may be relevant to study the production of a 
non-perturbative quark-gluon plasma at RHIC and LHC. 
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I. INTRODUCTION 



An exact non-perturbative result for electron-positron pair production from a constant 
electric field was obtained by Schwinger in 1951 by using proper time method In QCD 
this result depends on two independent casimir/gauge invariants C\ = [E a E a ] and C 2 = 
[d ahc E a E h E c f with color indices a, b, c=l,2,...8 in SU(3) |2j. Recently we have extended this 
calculation to arbitrary time dependent electric field E(t) in QED 3j and to arbitrary time 
dependent chromo-electric field E a (t) in QCD for quark-antiquark case 4|. This result relies 
crucially on the validity of the shift conjecture, which has not yet been established. In this 
paper we will extend this calculation to study non-perturbative gluon pair production from 
arbitrary E a (t). Unlike constant field E a case we encounter a non- vanishing single gluon 
(or tadpole) term in the presence of time dependent chromo-electric field E a {t). However, we 
show that the non-perturbative tadpole effective action d t^p T contains 5^ (pr) distribution 
and hence, for any non-vanishing transverse momentum, it does not contribute to non- 
perturbative gluon pair production rate via Schwinger mechanism. This result may 
be relevant to study the production of a non-perturbative quark-gluon plasma at RHIC and 
LHC Q, 5, 6, 7|. 

We obtain the following exact non-perturbative result for the probability of gluon (pair) 
production per unit time, per unit volume and per unit transverse momentum from an 
arbitrary time dependent chromo-electric field E a (t) with arbitrary color index a=l,2,...8 in 
SU(3): 

= i g 11111 + e '^ ] - (1) 

In the above equation 

A?(t) = ^[1 - cos0(t)]; A*3(t) = ^[1 + cos(^ ± e(t))}; cos*d{t) = -1 + 6C 2 (t)/Cf(t) 



(2) 



where 

C x {t) = [E a {t)E a {i)l and C 2 (t) = [d abc E a {t)E\t)E c {t)f (3) 

are two independent time-dependent casimir/gauge invariants in SU(3). 
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This result has the remarkable feature that it is independent of all the time derivatives 
d fL^ and has the same functional form as the constant chromo-electric field E a result 
with: E a -> E a (t). 

We will present a derivation of eq. (CQ) in this paper. 



II. SCHWINGER MECHANISM IN QCD IN THE PRESENCE OF ARBITRARY 
TIME-DEPENDENT CHROMO-ELECTRIC FIELD 



In the background field method of QCD [8|, |9[ the gauge field is the sum of classical 
chromo- field A a ^ and the quantum gluon field Q a . The non-abelian field tensor becomes 

F; U {A + Q] = d,(A a u + Ql) - dMl + QD + 9f ab \Al + Ql){Al + Ql). (4) 

The gauge field Lagrangian density is 



4* = -^ U [A + Q}F^[A + Q] -— [D,>[A\Cr? (5) 

where the second term in the right hand side is the gauge fixing term jfcj, . The covariant 
derivative is given by 

Df[A] = 5 a % + gf*°A°. (6) 

Keeping terms up to quadratic in Q field (for gluon pair production) and using Feynman-t' 
Hooft gauge (at=l) we find from eq. (jSJ) 

J d'xC = \jd A x \-{D,[A]Ql)F^[A] + Q^M^ U [A}Q» b ] 

= \jd 4 x [(D,[A]F^ a [A})Q a v + Q^Mf v [A]Q vb ] (7) 



where 



Mf v [A] = g, u [D 2 {A)r b - 2 9 r bc F c \A] (8) 



with g^ = (1, -1, -1, -1). 

The vacuum-to-vacuum transition amplitude for gluon (we will consider the ghost later 
in the derivation) in the presence of classical chromo-field A ^ is given by 

Z[A] fldQ}^ I d4x iQ^ M fii A ]Q" b MDAA]F^ a {A])Q<z] 
< °'° > ~~ Z[0] ~ JldQle* f d4x Q" aM $WQ vb ' ^ 
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We choose the arbitrary time- dependent chromo-electric field E a (t) to be along the z— axis 
(the beam direction) and work in the choice A% = [4| so that 

A a ,(x) = -6, E a (t)z. (10) 

The color indices a— 1,2, ...8 are arbitrary. From eq. ffTUl) we find 

D,[A\F^[A\ = 81 (11) 

Hence the single gluon term (D fl [A]F tl ' ya [A])Q® in eq. ([7j) was absent for constant chromo- 
electric field E a case in [3[ . However, in the presence of time-dependent chromo-electric field 
E a (t) this single gluon term in eq. ([7j) is not zero. Hence the path integration in eq. 
becomes more complicated due to the presence of this single gluon term. 

To evaluate the path integration in eq. Q we proceed as follows. We write 

M;l(x,x') = 5^(x-x')M; b L ,(x') (12) 

where we denote M°^(x) = M"^[A](x) which is given by eq. ([8]). The Green's function 



ib 



G ab (x,x')(= [M^Y^ix, x')) is given by 



d 4 x"M> lX ' ac (x, x")Gt(x", x') = 5 ab g^5 (4 \x - x'). (13) 

We change the variable 

<2» = Q'l( x ) ~ \ J d 4 x'Gf v {x,x')D x {x')F Xvh {x') (14) 

where we denote D^(x) = D^ b [A](x) and F®(x) = F^ v [A]{x). Under this change of variable 
[dQ] = [dQ'\. Using eqs. ATI]) . (fl3l and (P|l we find 

[{D^[A]F^ a [A))Q a v + Q» a M^[A}Q» b ] 

= -~ J d 4 x'D,{x)F» Xa {x)Gf v {x,x')D a {x')F™ h {x') + Q» a Mf v \A\Q' vb (15) 

Hence we find from eq. ([9]) 
a Z[A] 



< >' 



Z[0] 



exp[-'-Jd 4 x J d 4 x'D tl (x)F^ a {x)Gf v {x,x , )D a (x')F^ b (x')} 



x I = g-iStad x et = -iS tad asm (16 \ 
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where 



Stad = \ Jd A xJ d 4 x'D^x)F^ a (x)G 1/ x ab (x,x')D' T (x')F crvb (x') (17) 



is the tadpole (or single gluon) effective action and 



S {1) = -iln[ 



Bet- 1/2 M^[A] 



Det-V2M$[0] 
is the one loop (or gluon pair) effective action. 



(18) 



A. Tadpole (or Single Gluon) Effective Action in Arbitrary E a (t) 

For an operator M we will use Schwinger's notation [1] for the Green's function 
G(x,x') =< x\ — \x' >=< x\ / ds e~ sM \x' > . 
The Green's function G"^ u (x,x') for the operator M^ ab [A] becomes 

G» ab (x, x') =< x\ / ds e~ sM ™ W \ x > > . 
Using eq. ( fTOl) in ([H]) we find 



(19) 



(20) 



M™ b [A] = M a ; x [A]g^ = -5;[(5 ab p - tgf abc E c (t) Z y - 5 ab p 2 z - 5 ab p 2 T ] - 2gf abc E c (t)F»(2l) 
where 

r o 1 





l 1 

By using eq. ( 12T1) in ( !20l we find 

M ' Jo 
We write eq. (1231) in the Lorentz and color matrix notation as follows 



F v 



(22) 



> . (23) 



G^ fc (x,x')=y o d S [<x\ J ds e -(-(P0-9Mt)z) "+g4#+ 



(24) 



where 



A ab (t) = if abc E c (t) 



(25) 



Using eqs. fl2^|) and (iTTj) in (|T7|) we find the tadpole effective action 

^ = 1 f d'xd'x' Hds ^!M[< x | e -(-(P0-,A(t),) 2 + ^ + ^ + 2i 9 A(t)F) lx> >] 3a6^ 6 (i' 



2 J Jo dt 1 1 ^ dt 



1 f f 00 
- y d 2 :r T d 2 £y dt dz (it' dz' y ds 



dE a (t) i -s(-(pn-inA(t)z) 2 +v 2 ,+vl+ 2igA(t)F)\^/ ^ | _/ ^ u/ ^ 1 3a& ^ & 



dt L 1 1 1 1 T 1 1 J3 dt' v ; 

Inserting complete set of \px > states (by using J d 2 pT\pr >< Pt\ = 1) we find 

If , , , f°° dE a [t) 

Sta,& — - / d 2 xr d 2 x' T dt dz dt' dz' d 2 px / ds — : — [< xt\pt > 
2 J Jo dt 



dt' 



Using < q\p >= -j^e iqp we obt 



am 



Stad — —, — rrr / d 2 xx d 2 x' T dt dz dt' dz' d 2 px / ds 

2(2ir) 2 J Jo dt 

• T -|3afe^_ 

dt 



dE a (t) 



<t\ < r| e - s (-(P0-^ A W 2 ) 2 +Pz+PT+ 2igA{t)F)\ z i |./ -»x^-p T i3afe dE b (t') 

ii ii J3 r,, • v ; 



Integrating over x' T (by using J d 2 x' T e ix 'tVt = (2-n) 2 5^ 2 \pt)) we find 
dStzd 1 



- 5 (2) (p T ) e MT - pT y dt' y d^' y°°ds 



d 4 xd 2 pT 

dE a (t) ^ ^ 2 | e _ s( _ ( p _ i9 A (t)2) 2 + p2 +p 2 + 2i 9 A(t)F)| z / > | t / > pafc^ b ( t ') 



dt L " ~'~ '~ ' ' J3 dt 



(29) 



Since the above equation contains Dirac-delta function 5^ 2 \pt) we find (for any non- 
vanishing pt) 

,^! d = 0. (30) 

d 4 xd 2 p T K ' 

Hence the tadpole (or single gluon) effective action does not contribute to the exact result 
( d 4 xcPp T ) f° r ^ ne probability of gluon (pair) production per unit time per unit volume per 
unit transverse momentum from an arbitrary E a (t) via Schwinger mechanism. 

Now we evaluate the one-loop effective action for gluon and ghost in the presence of 
arbitrary E a (t) in the following. 
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B. One Loop (or Gluon Pair) Effective Action in Arbitrary E a {t) 

The one loop (or gluon pair) effective action eq. (fT8|) can be written as 

T)p\~ l l 2 M ab \A\ 7 
S W = -iW 1/9 „;7 r ] = -Ti[\nM^ b [A] - InilrtO]] 

= -Tr /°° ^£[ e ^(^r i 'M+«) - ^W^M-He)! ^ 31 ^ 
2 Jo s 

The trace Tr is given by 

TrO = tr Lorentz tr co i or / d 4 x < > . (32) 



Using eq. f[2"Tj) in f)3ip and writing in matrix notations we find 

i f ,a f°° ds 

'o s 



S W = -tr Lorentz tr color / d*X f°° —[< x \ e -^-9Ht)z?-vl-vl- 2i 9 K{t)F-ie) _ e - lS &-*)\ x >U 

2 J Jo s 

(33) 



where the Lorentz matrix and the color matrix A ab matrices are given by eqs. (|2"2"|) and 
( |25l) respectively. 

Using the eigen values 

-^eigenvalues = (^1; ^2; ^3) ^4) = (1; 0, —1, 0). (34) 

we perform the Lorentz trace and find 

= - tr co i or [J — J dt < t\ J dx < x\ J dy < y\ J dz < z\ 

e -is({p ~gh{t)z) 2 -pl~p 2 T ~ 2iX t gA(t)-ie) _ e ~is(p 2 -ie) |^ > |y > | x > |£ >] afe _ (35) 

Inserting complete set of \pr > states (using / d 2 px \pr >< Pt\ = 1) we find from the above 
equation 

[| +0 °c/t < t| / + °°^ < ^le-^w-^W*) 3 -^-*^*)]^ >\t>- J dt J dz^\} ab 

(36) 

where we have used the normalization < g|p >= -j=e iqp . At this stage we use the shift 



theorem 10j and find 

<? ' /*oo rj q p r o /"+oo /*+oo 



< z + 'd\ e ^KHt)^-vi-2^)u z + * 4- >\t>- I dt I dz^]r b 



gk{t)dt gk{t)dt J J Atts 

(37) 

where the z integration must be performed from — oo to +00 for the shift theorem to be 
applicable. 

Note that a state vector \z + ^y^- > which contains derivative operator is not familiar 
in physics. However, the state vector \z + > contains the derivative j| not ^. Hence 
the state vector is defined in the z-space with — acting as a c-number shift in z. To see how 
one operates with such state vector we find 

< z + -J7ST>* > /(*) = 4=e i(z+ ^)^ )p VW = -±=e^e~^f{t). (38) 
a^tj ar y An \J In 

Inserting complete sets of \p z > states (using / dp z \p z >< p z \ — 1) in eq. fl3Tl) we find 



S,,> = ^Sfj/fe/^'PM-/*/* (39) 



where 



= 7^ytr color [ £^°° dt<t\J dp z J dp' z £j dz e^e-imfiv* 

K p^ e is[-g 2 A' 2 (t)z 2 +pi+2iX ta A(t)]^r > e jK[t)ftP^ e - iz v'z\^ >]aft_ (4Q) 

1 d , 

It can be seen that the exponential e bW) contains the derivative 4 which operates on 
< p z |e ls ' _92A2 *'*'' z2+ ^' - * 2 ' Xi9A ^l|p2 > hence we can not move e~ sMt) mP z ^ q Yigixt. We insert more 
complete sets of states to find 

Fi{s) = 7^y tr coior [J dt J dt' J dt" J dz' J dz" J dz J dp J dp' Q J dp'^ J dp" 

J dp z J dp' z < t\p Q > e iPz < poe"^^ Pz \p' >< p' \t' >< t'\ < p z \z >< z'\ 
e is[- g ^K\t)z^+pl+2iX l9 K{t)]\ z " >< z "\ p ' z > \ t " >< t "\pl >< j%\e*®&'W > e~ lzp ' z < Po\t >T b 

j r+co f f f r+co f f f f 

= j2ny tTcolor 7 dt J dt' J dz' J dz" J dz J dp J dp' Q J dpi J dp" 

J dp z J dp' z e itpo e ipz <p Q e~^W)& Pz \p' > e - u ' p °e iz ' p * 



It can be seen that all the expressions in the above equation are independent of t except 
e it(po-Po"). This can be seen as follows 

< Po\f(t)j t \p' >= J dt' J dt" J dj%' < Pol*' >< t'\f{t)\t" >< t"\rt" >< rt"\j t \p' > 



dt' j dt" J dp^e- it,po 5(t' - f)f(i!')e*W ip' 5(p%'-p' ) = ip' I dt' e^ '^ f(f) 

(42) 

which is independent of t and 4. Hence by using the cyclic property of trace we can take the 
matrix [< p^e^tJ ^ \ p '£ >] ab to the left. The t integration is now easy (/+~ dte it( - p °- p '^ = 
2tt5(pq — p'q)) which gives 

Fi{s) = ^ptW [jdt'J dz' J dz" j + ™dz J dp, J dp' J dpi J d Pz J dp' z e Upo e izp * 

< p »\ e lAv&*\p >< p oe -^^\p' > e - iz ' p >e- lt ' p '° < z '\ e H^{t')z^+p%+2iX l9 A^) ]][Z n > 
e it'p^ e iz"p' ze -izp'^ab^ 

As advertised earlier we must integrate over z from — oo to +00 for the shift theorem to be 
applicable [uij]. The matrix element < z '\ e is [-9 2 ^{t)^+pl+^Mgk{t)}^ z ii > j s independent of z 
variable (it depends on z' and z" variables). Hence we can perform the z integration easily 
by using /+~ dze iz ^~ p '^ = 2n5(p z - p' z ) to find 



itpo 



Fi(s) = ^p tr coior [J dt' J dz' J dz" J dpo J dp' Q J dp" Q J dp z e 
< p'l\e^m^\p >< p e-^m^\p' > e-^e- 1 ^' < ^\ e is[- 9 ^Ht'> 2 +P 2 x +2iX ig A(t')]^ > 



Using the completeness relation / dp \p >< p \ = 1 we obtain 

Fi(s) = ^~p tr coior [J dt' J dz' J dz" J dp' J dp z 

e -iz'p z < z ^ e is{-g 2 A 2 (t')z 2 +p 2 z +2i\igA(t')]^ z ii > e iz"p z jab _ 



Since < z'\e ls ^ 9 2 A 2 (t')z 2 +p 2 +2iXigA(t')]^ z n > - g j nc l e p en( i en t of z variable (it depends on z' and 
z" variables) we can integrate over p z . We find (by using / dp z e^ z "~ z '^ Pz = (2tt)5(z' — z")) 

= 7^ytr color [JdtJdpoJ dz' < z '\ e ^-9 2 A 2 {t)z 2 +p 2 z +2iX ig A{t)} \ z f > ]a6 i ( 4g ) 

The Lorentz force equation in color space (in the adjoint representation of SU(3)) is given 
by 5 ab dp^ = gT^ b F^ u dx u = igf abc F^ v dx u . When the chromo-electric field is along the z-axis 
(eq. ( flOl) ) this becomes S ab dpo = igf abc E°(t)dz = gh ab {t)dz. Using this in ( H6l) and using 
the eigen values of the color matrix A ab (t) [2] from eq. ([2]) we find 

Fi(s) = j^T,Jdt gAj(t) J dz J dz' < z '\ e is{-9 2 A%t)z^+pl+2iX ig A j {t)]\ z ' > _ (4?) 
\ ) j=i 



The above equation boils down to an usual harmonic oscillator, u 2 (t)z' 2 + p 2 ,, with the 
constant frequency ui replaced by time dependent frequency oj{t). The normalized wave 
function is given by 4| 

J dz'\ < z'\n t > | 2 = 1. (48) 

Inserting complete set of harmonic oscillator states (by using J2n \ n t >< n t\ = 1) in eq. ff4T|) 
we find 

F i( s ) = T^-yl^Yl J dt 9&j{t) I dz J dz' < z'\n t > 



< n t \z >= j^rY^Y, J dt 9^j{t) J dz J dz'\ < z'\n t > \ 2 e 



e -sgA j (t){2n+l)+2X l gA j (t) 



2 -sgA J (t)(2n+l)+2\ l gA j (t) 



n j=\- 



1 6 . . sg2\ l Aj(t) 

Y, dt dz gA j (t) n , (49) 



(2n)f^J J JW 2sinh(sgAj(t)) 

where we have used eq. (I48p . Using this expression of Fi(s) in eq. (1391) and summing over I 
(by using the eigen values of the Dirac matrix from eq. ( 134}) ) we find 

167r 3 ^io si J F [y JW sinh(s^Aj(t)) s J v ' 

C. One Loop Effective Action for Ghost in arbitrary E a (t) 

Now we discuss the ghost contributions. The ghost Lagrangian density due to the gauge 



fixing term in eq. (jSJ) is given by 



0,3 



£gh = x a *[D„[A]D»[A + Q]]° b X b = X a] K ah [A,Q] x h (51) 

where x a is the ghost field. The vacuum-to-vacuum transition amplitude for ghost (anti- 
ghost) in the presence of A a (x) is given by: 

< 010 > A = M = Mjx[^^ = Detif^] = ^ 

1 " Z[0] /[dxt][d x ] e < /**x ta * afc[o1 *' DetAT^O] ' 1 ' 



This gives 



5 S = ~ Zln[ Det^[0] ] = -^t^ ab H - ln^[0]] (53) 

where X a6 [A] = K ab [A,Q = 0] is given by eq. (j5ip . We can repeat the ghost calculation 
similar to the gluon except for the following changes. There is an over all factor \ (with 
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[1 + cosh(2g A j(t))] — > 2) from eq. (150)) because we do not have any Lorentz matrices in 
if ab [A] in eq. (15T)) . There is another over all factor (—2) from eq. fl50|) because of the ghost 
determinant (compare eqs. (JTSJ) and (1551) ). With these two changes we find from eq. (150)) 
the following expression for the ghost one-loop effective action 

gh 167r 3 ^io si 7 w L sinh(s^A 3 (t)) s J v 7 

D. Non-Perturbative Gluon Pair Production From Arbitrary E a (t) via Schwinger 
Mechanism 

Adding eqs. (150)) and (1531) for the effective action for gluon and ghost respectively we 
find 

S$> = 77^ E f - />* /^^^[^(t) - -]■ (55) 

g 16vr 3 £^7o s i J F [S JW sinh sqA. t s J v ; 



The imaginary part of the above effective action gives real gluon pair production. The real 
part of the above equation is infrared divergent as s — > oo. However we are interested in 
the imaginary part of the above effective action which is not infrared divergent as s — > oo. 
This can be easily checked by making s — > is. The s-contour integration can be done in the 



similar way as was done in , lllj . Using the series expansion 



- + 2-E4^ (56) 



sinhx x ^ x 7i 2 n 2 + x 



we perform the s-contour integration around the pole s = ^° ^ nc ^ 

1 3 OO / 1 \n+l ~ ~ nnp^, 

W = 2lmS^ = -^E / d A x / d 2 p T \gAj(t)\e~^ . (57) 

47r 3=1 n=i n J J 

Hence the probability of non-perturbative gluon (pair) production per unit time, per unit 
volume and per unit transverse momentum from an arbitrary time dependent chromo-electric 
field E a (t) with arbitrary color index a=l,2,...8 in SU(3) is given by 

dW 1 3 np T 

- ^ E + (58) 

which reproduces eq. (IT)) . The expressions for gauge invariant Aj(£)'s are given in eq. (T2)). 
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III. CONCLUSION 



To conclude we have studied Schwinger mechanism for gluon pair production in the pres- 
ence of an arbitrary time- dependent chromo-electric background field E a (t) with arbitrary 
color index a=l,2,...8 in SU(3). We have obtained an exact result for the probability of 
non-perturbative gluon (pair) production per unit time per unit volume per unit transverse 
momentum d 4^p PT from arbitrary E a {t) by directly evaluating the path integral. We have 
found that the tadpole (or single gluon) effective action does not contribute to the non- 
perturbative gluon pair production rate d 4^ PT ■ We have found that the exact result for 
non-perturbative gluon pair production is independent of all the time derivatives d 
where n = 1,2, ...oo and has the same functional dependence on two casimir invariants 
\E a (t)E a (t)} and [d abc E a (t)E b (t)E°(t)) 2 as the constant chromo-electric field E a result with 
the replacement :E a — > E a (t). This result relies crucially on the validity of the shift conjec- 
ture, which has not yet been established. This result may be relevant to study the production 
of a non-perturbative quark-gluon plasma at RHIC and LHC 
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